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Abstract 

We propose studying the mechanism of factorization in exclusive decays of the 
form S ^ L>(*)X by examining the differential decay rate as a function of the 
invariant mass of the light hadronic state X. If factorization works primarily 
due to the large Nc limit then its accuracy is not expected to decrease as the 
X invariant mass increases. However, if factorization is mostly a consequence 
of perturbative QCD then the corrections should grow with the X invariant 
mass. Combining data for hadronic tau decays and semileptonic B decays 
allows tests of factorization to be made for a variety of final states. We 
discuss the examples oi B ^ D*-k^it^-k^'it^ and B D*ujit^. The mode 
B D*cvTr~ will allow a precision study of the dependence of the corrections 
to factorization on the invariant mass of the light hadronic state. 
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An understanding of nonleptonic B decays is important for the study of CP violation 
at B factories. However, there is at present no systematic, model-independent treatment 
of all such decays. In certain kinematic situations, factorization, according to which the 
matrix element of four-quark operators can be written as the product of pairs of matrix 
elements of two-quark operators, can be justified using perturbative QCD An example 
is B ^ D^*^X, where X denotes a hadronic state with low invariant mass, and the D*-*-* 
inherits the light degrees of freedom from the B meson. The large Nc limit of QCD |^ 
provides a different explanation of factorization, which is independent of the kinematics of 
the final state. Factorization has been shown to be consistent with experiment in certain two- 
body decays, such as B ^ D^*^n and D^*^p. Currently, this represents a test of factorization 
at the 10-20% level [Q. Factorization has also been observed to hold within the (presently 
sizable) errors in decays of the type B D^*^D^*\ where no perturbative QCD justification 
has been presented 0. 

In this paper we examine how the decays B D^*^X, where X is a final state containing 
two or more light hadrons, allow us to study the corrections to factorization^, and help 
determine the role that perturbative QCD plays. In such decays, factorization can be studied 
as a function of mx, the invariant mass of the final hadronic state produced by the light 
quark current. The perturbative QCD arguments for factorization depend on the light 
quarks being produced in an almost collinear state (for this reason, it has also been argued 
that factorization should hold for B D'^*^ +jet @]). Corrections are therefore expected to 
grow with mx/Ex, where Ex is the energy of the hadronic final state in the B rest frame, 
since this quantity characterizes the deviation of the hadronic state X from the light cone. 
However, if factorization works primarily due to the large Nc limit then its accuracy is not 
expected to decrease as mx increases. The behavior of factorization violations with mx 
can therefore distinguish the perturbative QCD explanation of factorization from that of 
the large Nc limit. Naively, the 1/N^ terms that violate factorization are 0(10%) in the 
amplitude and 0(20%) in the rate. When factorization is tested with a greater precision 
than this, it may be possible to see corrections to factorization which grow with mx, which 
would provide evidence that perturbative QCD plays an important role in factorization. 

Nonleptonic B meson decays arise dominantly due to the weak Hamiltonian 
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ci(mfe) (c7^Pi h) {qi^'PL u) + C2(mfe) {q^^Pl h) {cj^'Pl u) , (1) 



where Pl = (1 — 75)72 and q = d,s. The Hamiltonian in Eq. ([l|) is renormalized at the scale 
/i = mb- The coefficients q are given by Ci^2 = (c+ ± c_)/2, where 
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(2) 



in the leading logarithmic approximation. At next-to- leading order, ci(mb) = 1.13 and 
C2(^b) = —0.29 for as{Mz) = 0.118. Under the factorization hypothesis, the hadronic 
matrix element of H"^^ can be written as 



^These decays have also been discussed as tests of factorization by Reader and Isgur 
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(3) 



where the eUipses denote factorization violating terms. The first matrix element can be 
determined from semileptonic B — > D^*^ii' decays, while the second can be determined from 
the hadronic tau decay r Xur- 

The differential rate for tau decay to a given hadronic final state X may be written as 



drrix 327r2 m-^ 



"'^l {ml-mlf{ml + 2ml)vx{mj,), (4) 



where is the invariant mass-squared of the hadronic final state X, and vximj^) is related 
to the transverse part of the vacuum polarization amplitude. This holds in general even if 
X is a nonresonant multibody final state. The only assumption is that the weak hadronic 
current is conserved. For the vector part of the current, corrections to this assumption are 
suppressed by m^— m^, while for the axial part they are suppressed by rriu+md. Perturbative 
electroweak corrections have been neglected in Eq. (Q). 

The CLEO collaboration has measured fx(^x) for two-, three-, and four-pion final states 
QJ^. The two-pion final states are quasi- two-body, being dominated by the p resonance. 
Since the p is narrow, dr(i? — > D^*^ rnr) / drnj^ is a steeply falling function of in the region 
we are interested in, and so is less useful for our purposes than the three- and four-pion final 
states. The three-pion final states receive a large ai contribution, however, the ai is rather 
broad and there is evidence for nonresonant contributions in r decay to three pions. The 
four-pion final states are not dominated by any one resonance. 

In terms of fx(^x), the factorization prediction for the decay B D'^*'>X is 

^ " --?>Ti\c^{m,) + ^^^] vx{m\){l + 5NF). (5) 



X 



where (5nf denotes the nonfactorizable contributions, and in the denominator of the left- 
hand side mx denotes the dilepton invariant mass. In the large N^. limit d^p is of order 
whereas in the perturbative QCD approach it is expected to contain both perturbative 
corrections and nonperturbative ones suppressed by powers of mx/Ex, where Ex = {jn?B + 
m\ — m^(«))/(2m5) is the energy of the state X in the B rest frame. The ratio mx/Ex 
changes from 0.43 at mx = 1 GeV to 0.70 at mx = mr, so if there were order mx/Ex 
corrections to factorization then the accuracy of the predictions should change significantly 
over the accessible range of mx- 

In Fig. |l| we plot the distribution in B ^ D*tt^ 71^71^71^ decay, normalized to 
the B — > D*ih' rate, as predicted by Eq. (|^). We used the CLEO fit to the shape of 
dr[B D*iu)/dm'j^ [jlO| and the r decay data 0. This decay mode is interesting for test- 



ing factorization because it is not dominated by a single narrow resonance, and its branching 



fraction is large, almost 2% |TTI|. However, this mode has a significant disadvantage for test- 
ing factorization because of a potentially sizeable background from processes where one or 
more of the pions are emitted from the (cb) current. In this case, the (cb) current creates a 
charm state C of the type D*-|-n7r(l<n<3), and the {du) current creates the remaining 
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FIG. 1. dT{B — > D* t:^ TT^ -n^ t:^) / dm\ , where mx is the 1:^1:^1:^11^ invariant mass, normaUzed 
to the semileptonic width T{B D*£i/). The triangles are the B decay data [11,1^, and the 
squares are our predictions using the r data M. There is an additional 9% uncertainty in the B 



decay data from the overall normalization [12|. 



4 — n pions in the final state. In the formal limit when ruh — rUc is large this background is 
parametrically suppressed, since the light hadrons from the C decay combined with those 
from the {du) current will have large invariant mass over most of the phase space. This can, 
however, be a sizable background for the physical rrih and rric masses. 

Because of the problem of backgrounds of this type, a more promising way to test the 
m'j^ dependence of factorization is to focus on the subset of these decays which proceeds 
via B — >■ D*uJir~. In this case, the backgrounds are easier to constrain, because the (du) 
current must create the vr" in the final state. Therefore, the (cb) current would have to 
create charmed states, C(^d*^'), which can decay into D*u!. This background can be related 
using factorization (in the most reliable case of the current only making one pion) to the 



B -^C, 



decay rate. Assuming that C(£).(^) is a narrow resonance. 
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dw 
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where w = vb ■ vc- Assuming that the C(^£)*uj) niass is its minimal value m™™ = m^* + m^^ 
and that semileptonic B — > C^o'u))^^ rate depends on w as ^/w"^ — 1, we find that V{B ^ 
C{D*uj)T^) — 0.1 X V{B C{D-'u})^^)- Therefore, this background can be constrained by 
measuring the B D*ujiu decay rate. The measured B D*ujt:^ branching fraction is 
about 0.3% [|n|]. Thus, for example, a bound of B{B Ci^o'u))^^) < 0.3% would imply that 
the feed-down from higher mass Ci^r,*ui) states to the D*ujir^ final state is less than 10%. 
This is a very conservative estimate, since the Ci^n-'u)) branching fraction to D*uj is probably 
significantly less than unity. Moreover, the B — >■ C(^D*^)iiJ branching fraction is probably 
much smaller than 0.3%; this is certainly the case in most models |13]. Therefore, we expect 
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FIG. 2. dr(i? — > L'*a;7r)/dmj^, where is the wtt" invariant mass, normahzed to the semilep- 
tonic width Y{B — > D*ii>). The triangles are the B decay data |11|, and the squares are our 
predictions using the T data g. There is an additional 9% uncertainty in the B decay data from 



the overall normalization [12|. 



that the B D*ujt:~ decay will only receive a very small background from higher mass 
charm states. 

In Fig. 1^ we plot the m\ distribution m B ^ 0*0011^ decay, normalized to the B —>■ D*li> 
rate, as predicted by Eq. (^) and the r decay data 0. This decay receives a large contribution 
from a p' resonance with mass around 1.4—1.5 GeV P,pT|. However, the p' is very broad, with 
a width estimated to be around 400 — 500 MeV. Since we would like to make a precision test 
of factorization and learn about the size of corrections suppressed by powers of mx / Ex and 
Q!s, one needs tests that do not make assumptions based on resonance models or heavy quark 
symmetry [|14|, which may have corrections of similar size. This is best achieved by looking 
at the differential mass distribution, which should also teach us about how factorization 
breaks down (if it does at all) at higher masses. 

Controlling this background in the B D*7r+7r^7r^7r° decay is more complicated than 
for B — i> D*uJ7r~ , and may make it difficult to make a precise study of factorization in this 
channel. The contribution when C contains two or three pions is expected to be small and 
can be constrained by bounding the semileptonic rates as in the case of B D*ujtt^. The 
contribution when one pion comes from the {cb) current and three pions from the {du) current 
may be estimated from the measured f37r(m^) and the B — > D*iTiu branching fraction. 
Using the ALEPH data |jT5| for f37r(m|^) measured in r decay and B{B — >• Z)*7r£z/) ~ 1%, it 
appears that this background may be significant. The size of this background can be tested 
experimentally by measuring the B —>■ Z}*°7r+7r+7r^7r^ decay rate, since this decay only 
receives contributions from higher mass charm states which decay into a D*^ plus one or three 
pions. Unless the B — > Z}*°7r"'"7r"'"7r~7r~ rate is much smaller than the B D*ti^tt~ti~tt^ 
rate, it will be difficult to make a precise test of factorization in B ^ 0*71^71^71^71^ . 

Data on e~^e~ — >• hadrons can also be used to determine vx{Tnx) and predict the mass 
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spectrum in B ^ D^*^X decay |T^. Currently the CMD-2 data on e+e" An are available 
up to = (1.4 GeV)^ measurements of this cross section at higher energies would 
allow testing factorization at values of > m^. 

In summary, for a precision study of corrections to factorization in i? ^ D^*^X decays 
(where X contains only light hadrons), it is best not to make assumptions based on heavy 
quark symmetry or resonant decomposition (especially for broad states). Rather, one should 
use the differential mj^ distribution with a fixed hadron content and compare with the 
predicted rate based on factorization using the spectral function extracted from r decay. 
Over the accessible kinematic range the parameter rrix/Ex, which controls the accuracy of 
the predictions if the primary reason for factorization is perturbative QCD, changes from 
about 0.4 to 0.7, making mx-dependent violations of factorization potentially observable. In 
this paper we concentrated on final states containing four pions, but this approach is general 
and can be applied to test factorization in decays to any light hadronic final state where 
backgrounds from the (cb) current producing some of the light hadrons can be controlled. 
Studying the available data in B ^ D*7i^7i^7i^7i^ and B D*ujtt~ decays, we found no 
evidence of factorization becoming a worse approximation at higher values of mx- It will 
be very interesting to make a precision study, which should be possible at the B factories. 
The decay mode B D*ujti~ is particularly well suited for this purpose. 
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